Step 1: Generate Stress/No-Stress Labels: We replace the HR response variable with the binary variable Stress/No-Stress. For that, we generate Stress/No-Stress labels for each 5 min interval based on the HR values and the psycho-physiological definition of stress as instantaneous hyperarousal. If normalized HR does not exhibit severe deviations from normality, then anything to the right of mean could be considered hyperarousal (and thus Stress). The remaining of the distribution could be considered as No-Stress.

Step 2: Regress on Stress/No-Stress: Using the Stress/No-Stress binary variable as the response variable, we run a model that provides explanations for the sources of participant stress.


1 Age and BMI Descriptive Statistics

2 Physical Acitvity (PA) Per Participant & Day

3 Frequencies of Sedentary Activities Per Day and Participant

## , ,  = CD
## 
##       
##        Work NonWork Driving
##   T001  133       0      10
##   T002  131      41       0
##   T003  137       5      10
##   T004  147       3      12
##   T005  109      17      16
##   T006   74      42       0
##   T007  128       1       5
##   T008   61      27      20
##   T009   59      10       0
##   T010  100       0      11
##   T011  100      11      10
##   T012   68      87      19
##   T014   99      30       0
##   T015   87      26       6
##   T016  105      42       4
##   T017  131       0       0
##   T018   74      31       0
##   T019   66      18      50
##   T020  136      14      31
##   T021   38      82      40
##   T022  149      51       0
##   T023  159       1      15
##   T024  119      11       0
##   T025  127      13      35
## 
## , ,  = WD
## 
##       
##        Work NonWork Driving
##   T001  162      35       0
##   T002   88     145      23
##   T003  109      34     150
##   T004  344      16       8
##   T005  148     119      69
##   T006  194      54      38
##   T007  187      29      73
##   T008  195      16      18
##   T009  254       8       0
##   T010  164      64      16
##   T011   82      18      43
##   T012   75      75      10
##   T014  145     122      13
##   T015  211      50      19
##   T016  196      55      15
##   T017   96      42      20
##   T018  184      57      15
##   T019  106     184      20
##   T020  232      25      20
##   T021  103      85      63
##   T022  169     143       0
##   T023  199      57      52
##   T024  191      29      38
##   T025  223      67      86
## 
## , ,  = ND
## 
##       
##        Work NonWork Driving
##   T001    0     117       4
##   T002    0     122       0
##   T003    0     109      10
##   T004    0     160       8
##   T005    0     103      21
##   T006    0      76      15
##   T007    0     117       0
##   T008    0      87       6
##   T009    0      80       0
##   T010    0      98       0
##   T011    0     118      19
##   T012    0      98       7
##   T014    0     111      43
##   T015    0     161      13
##   T016    0     152       0
##   T017    0      52       4
##   T018    0     117       5
##   T019    0     168       6
##   T020    0      91      23
##   T021    0     109      29
##   T022    0      98       0
##   T023    0     112      15
##   T024    0      95       0
##   T025    0     141       0

4 Distrubution of NHR in Sedentary Activities

5 Viz of Stress Distribution and QQ Plots for T025

6 Viz of Stress Distribution per Participant and Day

7 PA_Percent Plot

8 Stress Percent Plots

9 Big 5 Box Plots

10 Cadence PDF

11 Profession and Gender Plot

12 NASA Mental Demand Plot

## 
##  Descriptive statistics by group 
## group: CD
##             vars  n  mean    sd median trimmed   mad min max range  skew
## Participant    1 24 12.50  7.07   12.5    12.5  8.90   1  24    23  0.00
## Day            2 24  1.00  0.00    1.0     1.0  0.00   1   1     0   NaN
## Date           3 24 28.88 20.24   27.5    28.1 25.95   2  65    63  0.21
## N_MD           4 24 10.67  5.74   11.5    10.8  6.67   1  19    18 -0.20
## Profession     5 24  1.42  0.50    1.0     1.4  0.00   1   2     1  0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.34 4.13
## N_MD           -1.35 1.17
## Profession     -1.98 0.10
## ------------------------------------------------------------ 
## group: WD
##             vars  n  mean    sd median trimmed   mad min max range skew
## Participant    1 45 12.47  7.21     13   12.46  8.90   1  24    23 0.01
## Day            2 45  2.00  0.00      2    2.00  0.00   2   2     0  NaN
## Date           3 45 30.51 21.57     34   29.89 31.13   1  69    68 0.13
## N_MD           4 45  8.00  4.40      8    7.68  5.93   2  17    15 0.52
## Profession     5 45  1.44  0.50      1    1.43  0.00   1   2     1 0.22
##             kurtosis   se
## Participant    -1.37 1.07
## Day              NaN 0.00
## Date           -1.39 3.21
## N_MD           -0.92 0.66
## Profession     -2.00 0.07
## ------------------------------------------------------------ 
## group: ND
##             vars  n  mean    sd median trimmed   mad min max range skew
## Participant    1 24 12.50  7.07   12.5   12.50  8.90   1  24    23 0.00
## Day            2 24  3.00  0.00    3.0    3.00  0.00   3   3     0  NaN
## Date           3 24 32.08 21.15   31.5   31.35 27.43   4  68    64 0.18
## N_MD           4 24  4.88  4.77    3.0    4.15  1.48   1  16    15 1.39
## Profession     5 24  1.42  0.50    1.0    1.40  0.00   1   2     1 0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.36 4.32
## N_MD            0.51 0.97
## Profession     -1.98 0.10

13 NASA Physical Demand Plot

## 
##  Descriptive statistics by group 
## group: CD
##             vars  n  mean    sd median trimmed   mad min max range skew
## Participant    1 24 12.50  7.07   12.5   12.50  8.90   1  24    23 0.00
## Day            2 24  1.00  0.00    1.0    1.00  0.00   1   1     0  NaN
## Date           3 24 28.88 20.24   27.5   28.15 25.95   2  65    63 0.21
## N_PD           4 24  7.46  5.34    5.0    6.95  3.71   2  18    16 0.78
## Profession     5 24  1.42  0.50    1.0    1.40  0.00   1   2     1 0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.34 4.13
## N_PD           -0.88 1.09
## Profession     -1.98 0.10
## ------------------------------------------------------------ 
## group: WD
##             vars  n  mean    sd median trimmed   mad min max range skew
## Participant    1 45 12.47  7.21     13   12.46  8.90   1  24    23 0.01
## Day            2 45  2.00  0.00      2    2.00  0.00   2   2     0  NaN
## Date           3 45 30.51 21.57     34   29.89 31.13   1  69    68 0.13
## N_PD           4 45  5.44  3.70      4    4.92  2.97   1  16    15 1.28
## Profession     5 45  1.44  0.50      1    1.43  0.00   1   2     1 0.22
##             kurtosis   se
## Participant    -1.37 1.07
## Day              NaN 0.00
## Date           -1.39 3.21
## N_PD            1.19 0.55
## Profession     -2.00 0.07
## ------------------------------------------------------------ 
## group: ND
##             vars  n  mean    sd median trimmed   mad min max range skew
## Participant    1 24 12.50  7.07   12.5    12.5  8.90   1  24    23 0.00
## Day            2 24  3.00  0.00    3.0     3.0  0.00   3   3     0  NaN
## Date           3 24 32.08 21.15   31.5    31.4 27.43   4  68    64 0.18
## N_PD           4 24  4.67  3.25    4.0     4.3  2.22   1  15    14 1.37
## Profession     5 24  1.42  0.50    1.0     1.4  0.00   1   2     1 0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.36 4.32
## N_PD            1.86 0.66
## Profession     -1.98 0.10

14 NASA Temporal Demand Plot

15 NASA Performance Plot

## 
##  Descriptive statistics by group 
## group: CD
##             vars  n  mean    sd median trimmed   mad min max range  skew
## Participant    1 24 12.50  7.07   12.5   12.50  8.90   1  24    23  0.00
## Day            2 24  1.00  0.00    1.0    1.00  0.00   1   1     0   NaN
## Date           3 24 28.88 20.24   27.5   28.15 25.95   2  65    63  0.21
## N_PD           4 24  7.46  5.34    5.0    6.95  3.71   2  18    16  0.78
## N_P            5 24 14.50  6.44   17.5   15.30  2.22   1  20    19 -1.17
## Profession     6 24  1.42  0.50    1.0    1.40  0.00   1   2     1  0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.34 4.13
## N_PD           -0.88 1.09
## N_P            -0.15 1.31
## Profession     -1.98 0.10
## ------------------------------------------------------------ 
## group: WD
##             vars  n  mean    sd median trimmed   mad min max range  skew
## Participant    1 45 12.47  7.21     13   12.46  8.90   1  24    23  0.01
## Day            2 45  2.00  0.00      2    2.00  0.00   2   2     0   NaN
## Date           3 45 30.51 21.57     34   29.89 31.13   1  69    68  0.13
## N_PD           4 45  5.44  3.70      4    4.92  2.97   1  16    15  1.28
## N_P            5 45 15.58  3.68     16   15.86  4.45   6  20    14 -0.53
## Profession     6 45  1.44  0.50      1    1.43  0.00   1   2     1  0.22
##             kurtosis   se
## Participant    -1.37 1.07
## Day              NaN 0.00
## Date           -1.39 3.21
## N_PD            1.19 0.55
## N_P            -0.55 0.55
## Profession     -2.00 0.07
## ------------------------------------------------------------ 
## group: ND
##             vars  n  mean    sd median trimmed   mad min max range  skew
## Participant    1 24 12.50  7.07   12.5    12.5  8.90   1  24    23  0.00
## Day            2 24  3.00  0.00    3.0     3.0  0.00   3   3     0   NaN
## Date           3 24 32.08 21.15   31.5    31.4 27.43   4  68    64  0.18
## N_PD           4 24  4.67  3.25    4.0     4.3  2.22   1  15    14  1.37
## N_P            5 24 12.83  4.91   12.5    13.2  4.45   2  20    18 -0.38
## Profession     6 24  1.42  0.50    1.0     1.4  0.00   1   2     1  0.32
##             kurtosis   se
## Participant    -1.35 1.44
## Day              NaN 0.00
## Date           -1.36 4.32
## N_PD            1.86 0.66
## N_P            -0.50 1.00
## Profession     -1.98 0.10

16 NASA Effort Plot

17 NASA Frustration Plot

18 Sleep Duration PDF

19 mCadence Plot

20 mCadence Box Plot

21 EDA Plots

22 Correlation Plot

23 Full Logistic Regression Mixed Effects Model (Stress from NHR)

## Generalized linear mixed model fit by maximum likelihood (Laplace
##   Approximation) [glmerMod]
##  Family: binomial  ( logit )
## Formula: NHR_Stress ~ Cadence + Day * Period + Activity4 + Day * PA_Percent *  
##     PA_Intensity + Sleep_Time + Profession + N_MD + N_PD + N_P +  
##     Agreeableness + Conscientiousness + Extraversion + Neuroticism +  
##     Openness + (1 | Participant)
##    Data: df_sedentMultimodal
## 
##      AIC      BIC   logLik deviance df.resid 
##    12575    12791    -6259    12517    12679 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -13.300  -0.608  -0.249   0.633   6.944 
## 
## Random effects:
##  Groups      Name        Variance Std.Dev.
##  Participant (Intercept) 1.17     1.08    
## Number of obs: 12708, groups:  Participant, 24
## 
## Fixed effects:
##                               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)                   -7.03651    2.59726   -2.71   0.0067 ** 
## Cadence                        0.13492    0.00402   33.58  < 2e-16 ***
## DayWD                          4.35547    0.72953    5.97  2.4e-09 ***
## DayND                          8.47083    0.65730   12.89  < 2e-16 ***
## PeriodAfternoon                0.54805    0.09449    5.80  6.6e-09 ***
## Activity4NonWork              -0.78149    0.06960  -11.23  < 2e-16 ***
## Activity4Driving              -0.06852    0.07928   -0.86   0.3874    
## PA_Percent                     0.64288    0.04630   13.89  < 2e-16 ***
## PA_Intensity                   0.18426    0.01563   11.79  < 2e-16 ***
## Sleep_Time                    -0.06051    0.02500   -2.42   0.0155 *  
## ProfessionSP                   0.73503    0.50572    1.45   0.1461    
## N_MD                           0.09079    0.00903   10.05  < 2e-16 ***
## N_PD                           0.08821    0.00921    9.58  < 2e-16 ***
## N_P                            0.01460    0.00666    2.19   0.0284 *  
## Agreeableness                  0.11423    0.14248    0.80   0.4227    
## Conscientiousness              0.06547    0.15059    0.43   0.6637    
## Extraversion                  -0.30581    0.14812   -2.06   0.0390 *  
## Neuroticism                   -0.10907    0.16912   -0.64   0.5190    
## Openness                      -0.19145    0.15353   -1.25   0.2124    
## DayWD:PeriodAfternoon         -0.36880    0.11600   -3.18   0.0015 ** 
## DayND:PeriodAfternoon         -0.67042    0.14413   -4.65  3.3e-06 ***
## DayWD:PA_Percent              -0.40581    0.05629   -7.21  5.6e-13 ***
## DayND:PA_Percent              -0.49641    0.05274   -9.41  < 2e-16 ***
## DayWD:PA_Intensity            -0.10852    0.01852   -5.86  4.7e-09 ***
## DayND:PA_Intensity            -0.20665    0.01695  -12.19  < 2e-16 ***
## PA_Percent:PA_Intensity       -0.01686    0.00115  -14.63  < 2e-16 ***
## DayWD:PA_Percent:PA_Intensity  0.00960    0.00138    6.97  3.1e-12 ***
## DayND:PA_Percent:PA_Intensity  0.01461    0.00129   11.32  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## optimizer (Nelder_Mead) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.625385 (tol = 0.002, component 1)
## Model is nearly unidentifiable: very large eigenvalue
##  - Rescale variables?
## Model is nearly unidentifiable: large eigenvalue ratio
##  - Rescale variables?
Observations 12708
Dependent variable NHR_Stress
Type Mixed effects generalized linear model
Family binomial
Link logit
AIC 12575.1707
BIC 12791.2204
Pseudo-R² (fixed effects) 0.3550
Pseudo-R² (total) 0.5246
Fixed Effects
Est. S.E. z val. p
(Intercept) -7.0365 2.5973 -2.7092 0.0067
Cadence 0.1349 0.0040 33.5816 0.0000
DayWD 4.3555 0.7295 5.9702 0.0000
DayND 8.4708 0.6573 12.8873 0.0000
PeriodAfternoon 0.5480 0.0945 5.8002 0.0000
Activity4NonWork -0.7815 0.0696 -11.2280 0.0000
Activity4Driving -0.0685 0.0793 -0.8643 0.3874
PA_Percent 0.6429 0.0463 13.8861 0.0000
PA_Intensity 0.1843 0.0156 11.7856 0.0000
Sleep_Time -0.0605 0.0250 -2.4205 0.0155
ProfessionSP 0.7350 0.5057 1.4534 0.1461
N_MD 0.0908 0.0090 10.0515 0.0000
N_PD 0.0882 0.0092 9.5797 0.0000
N_P 0.0146 0.0067 2.1922 0.0284
Agreeableness 0.1142 0.1425 0.8018 0.4227
Conscientiousness 0.0655 0.1506 0.4348 0.6637
Extraversion -0.3058 0.1481 -2.0645 0.0390
Neuroticism -0.1091 0.1691 -0.6449 0.5190
Openness -0.1914 0.1535 -1.2470 0.2124
DayWD:PeriodAfternoon -0.3688 0.1160 -3.1794 0.0015
DayND:PeriodAfternoon -0.6704 0.1441 -4.6516 0.0000
DayWD:PA_Percent -0.4058 0.0563 -7.2093 0.0000
DayND:PA_Percent -0.4964 0.0527 -9.4127 0.0000
DayWD:PA_Intensity -0.1085 0.0185 -5.8589 0.0000
DayND:PA_Intensity -0.2066 0.0170 -12.1891 0.0000
PA_Percent:PA_Intensity -0.0169 0.0012 -14.6320 0.0000
DayWD:PA_Percent:PA_Intensity 0.0096 0.0014 6.9715 0.0000
DayND:PA_Percent:PA_Intensity 0.0146 0.0013 11.3174 0.0000
Random Effects
Group Parameter Std. Dev.
Participant (Intercept) 1.0800
Grouping Variables
Group # groups ICC
Participant 24 0.2629
## $Participant
##      (Intercept)
## T001      0.6509
## T002     -0.7428
## T003      0.7691
## T004      0.0211
## T005     -0.8043
## T006     -0.9528
## T007      0.0074
## T008     -1.2787
## T009     -0.0232
## T010     -0.0912
## T011     -1.0734
## T012     -0.9816
## T014      1.1980
## T015      1.2066
## T016     -1.9259
## T017      2.1585
## T018      0.2552
## T019      1.8873
## T020      1.5528
## T021     -0.0905
## T022     -0.4335
## T023      0.0743
## T024     -1.7718
## T025      0.4348
## 
## with conditional variances for "Participant"

24 Remove Conscientiousness from the Model

## Generalized linear mixed model fit by maximum likelihood (Laplace
##   Approximation) [glmerMod]
##  Family: binomial  ( logit )
## Formula: NHR_Stress ~ Cadence + Day * Period + Activity4 + Day * PA_Percent *  
##     PA_Intensity + Sleep_Time + Profession + N_MD + N_PD + N_P +  
##     Agreeableness + Extraversion + Neuroticism + Openness + (1 |  
##     Participant)
##    Data: df_sedentMultimodal
## 
##      AIC      BIC   logLik deviance df.resid 
##    12573    12782    -6259    12517    12680 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -13.297  -0.608  -0.249   0.633   6.954 
## 
## Random effects:
##  Groups      Name        Variance Std.Dev.
##  Participant (Intercept) 1.18     1.09    
## Number of obs: 12708, groups:  Participant, 24
## 
## Fixed effects:
##                               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)                   -6.64120    2.55185   -2.60   0.0093 ** 
## Cadence                        0.13493    0.00402   33.58  < 2e-16 ***
## DayWD                          4.36338    0.73592    5.93  3.0e-09 ***
## DayND                          8.47634    0.66333   12.78  < 2e-16 ***
## PeriodAfternoon                0.54792    0.09454    5.80  6.8e-09 ***
## Activity4NonWork              -0.78169    0.06961  -11.23  < 2e-16 ***
## Activity4Driving              -0.06845    0.07928   -0.86   0.3879    
## PA_Percent                     0.64328    0.04662   13.80  < 2e-16 ***
## PA_Intensity                   0.18436    0.01576   11.70  < 2e-16 ***
## Sleep_Time                    -0.06068    0.02500   -2.43   0.0152 *  
## ProfessionSP                   0.77390    0.49995    1.55   0.1216    
## N_MD                           0.09068    0.00903   10.04  < 2e-16 ***
## N_PD                           0.08840    0.00920    9.61  < 2e-16 ***
## N_P                            0.01460    0.00666    2.19   0.0284 *  
## Agreeableness                  0.11861    0.14430    0.82   0.4111    
## Extraversion                  -0.29982    0.15085   -1.99   0.0469 *  
## Neuroticism                   -0.11865    0.17234   -0.69   0.4912    
## Openness                      -0.17977    0.15522   -1.16   0.2468    
## DayWD:PeriodAfternoon         -0.36856    0.11605   -3.18   0.0015 ** 
## DayND:PeriodAfternoon         -0.66979    0.14423   -4.64  3.4e-06 ***
## DayWD:PA_Percent              -0.40642    0.05665   -7.17  7.2e-13 ***
## DayND:PA_Percent              -0.49674    0.05305   -9.36  < 2e-16 ***
## DayWD:PA_Intensity            -0.10871    0.01867   -5.82  5.8e-09 ***
## DayND:PA_Intensity            -0.20678    0.01709  -12.10  < 2e-16 ***
## PA_Percent:PA_Intensity       -0.01687    0.00116  -14.54  < 2e-16 ***
## DayWD:PA_Percent:PA_Intensity  0.00961    0.00138    6.94  3.9e-12 ***
## DayND:PA_Percent:PA_Intensity  0.01462    0.00130   11.26  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## optimizer (Nelder_Mead) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.612856 (tol = 0.002, component 1)
## Model is nearly unidentifiable: very large eigenvalue
##  - Rescale variables?
## Model is nearly unidentifiable: large eigenvalue ratio
##  - Rescale variables?
Observations 12708
Dependent variable NHR_Stress
Type Mixed effects generalized linear model
Family binomial
Link logit
AIC 12573.3571
BIC 12781.9568
Pseudo-R² (fixed effects) 0.3532
Pseudo-R² (total) 0.5244
Fixed Effects
Est. S.E. z val. p
(Intercept) -6.6412 2.5519 -2.6025 0.0093
Cadence 0.1349 0.0040 33.5809 0.0000
DayWD 4.3634 0.7359 5.9291 0.0000
DayND 8.4763 0.6633 12.7785 0.0000
PeriodAfternoon 0.5479 0.0945 5.7956 0.0000
Activity4NonWork -0.7817 0.0696 -11.2291 0.0000
Activity4Driving -0.0685 0.0793 -0.8635 0.3879
PA_Percent 0.6433 0.0466 13.7987 0.0000
PA_Intensity 0.1844 0.0158 11.6957 0.0000
Sleep_Time -0.0607 0.0250 -2.4274 0.0152
ProfessionSP 0.7739 0.5000 1.5479 0.1216
N_MD 0.0907 0.0090 10.0424 0.0000
N_PD 0.0884 0.0092 9.6052 0.0000
N_P 0.0146 0.0067 2.1923 0.0284
Agreeableness 0.1186 0.1443 0.8219 0.4111
Extraversion -0.2998 0.1509 -1.9876 0.0469
Neuroticism -0.1186 0.1723 -0.6885 0.4912
Openness -0.1798 0.1552 -1.1582 0.2468
DayWD:PeriodAfternoon -0.3686 0.1161 -3.1758 0.0015
DayND:PeriodAfternoon -0.6698 0.1442 -4.6439 0.0000
DayWD:PA_Percent -0.4064 0.0566 -7.1747 0.0000
DayND:PA_Percent -0.4967 0.0531 -9.3630 0.0000
DayWD:PA_Intensity -0.1087 0.0187 -5.8215 0.0000
DayND:PA_Intensity -0.2068 0.0171 -12.0989 0.0000
PA_Percent:PA_Intensity -0.0169 0.0012 -14.5354 0.0000
DayWD:PA_Percent:PA_Intensity 0.0096 0.0014 6.9399 0.0000
DayND:PA_Percent:PA_Intensity 0.0146 0.0013 11.2566 0.0000
Random Effects
Group Parameter Std. Dev.
Participant (Intercept) 1.0900
Grouping Variables
Group # groups ICC
Participant 24 0.2647

25 Optimal and Scaled Model

## Generalized linear mixed model fit by maximum likelihood (Laplace
##   Approximation) [glmerMod]
##  Family: binomial  ( logit )
## Formula: 
## NHR_Stress ~ s_Cadence + Day * Period + Activity4 + Day * s_PA_Percent *  
##     s_PA_Intensity + s_Sleep_Time + Profession + s_N_MD + s_N_PD +  
##     s_N_P + s_Agreeableness + s_Extraversion + s_Neuroticism +  
##     s_Openness + (1 | Participant)
##    Data: df_sedentMultimodal
## 
##      AIC      BIC   logLik deviance df.resid 
##    12573    12782    -6259    12517    12680 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -13.296  -0.608  -0.249   0.633   6.955 
## 
## Random effects:
##  Groups      Name        Variance Std.Dev.
##  Participant (Intercept) 1.18     1.09    
## Number of obs: 12708, groups:  Participant, 24
## 
## Fixed effects:
##                                   Estimate Std. Error z value Pr(>|z|)    
## (Intercept)                        -0.7764     0.3167   -2.45   0.0142 *  
## s_Cadence                           0.9420     0.0281   33.58  < 2e-16 ***
## DayWD                              -0.2166     0.0988   -2.19   0.0283 *  
## DayND                               1.3800     0.1405    9.82  < 2e-16 ***
## PeriodAfternoon                     0.5482     0.0945    5.80  6.5e-09 ***
## Activity4NonWork                   -0.7817     0.0696  -11.23  < 2e-16 ***
## Activity4Driving                   -0.0684     0.0793   -0.86   0.3884    
## s_PA_Percent                        0.2893     0.0642    4.50  6.7e-06 ***
## s_PA_Intensity                     -0.0177     0.0657   -0.27   0.7877    
## s_Sleep_Time                       -0.0928     0.0382   -2.43   0.0152 *  
## ProfessionSP                        0.7743     0.4995    1.55   0.1211    
## s_N_MD                              0.4721     0.0470   10.04  < 2e-16 ***
## s_N_PD                              0.3489     0.0363    9.61  < 2e-16 ***
## s_N_P                               0.0698     0.0319    2.19   0.0285 *  
## s_Agreeableness                     0.2090     0.2528    0.83   0.4084    
## s_Extraversion                     -0.5591     0.2774   -2.02   0.0439 *  
## s_Neuroticism                      -0.1907     0.2723   -0.70   0.4837    
## s_Openness                         -0.3048     0.2589   -1.18   0.2392    
## DayWD:PeriodAfternoon              -0.3690     0.1160   -3.18   0.0015 ** 
## DayND:PeriodAfternoon              -0.6703     0.1441   -4.65  3.3e-06 ***
## DayWD:s_PA_Percent                 -0.4303     0.0867   -4.96  7.0e-07 ***
## DayND:s_PA_Percent                  0.1514     0.0881    1.72   0.0859 .  
## DayWD:s_PA_Intensity               -0.0224     0.0866   -0.26   0.7960    
## DayND:s_PA_Intensity               -0.3932     0.0902   -4.36  1.3e-05 ***
## s_PA_Percent:s_PA_Intensity        -0.9732     0.0682  -14.27  < 2e-16 ***
## DayWD:s_PA_Percent:s_PA_Intensity   0.5549     0.0812    6.83  8.2e-12 ***
## DayND:s_PA_Percent:s_PA_Intensity   0.8440     0.0765   11.03  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## optimizer (Nelder_Mead) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00221052 (tol = 0.002, component 1)
Observations 12708
Dependent variable NHR_Stress
Type Mixed effects generalized linear model
Family binomial
Link logit
AIC 12573.3561
BIC 12781.9558
Pseudo-R² (fixed effects) 0.3534
Pseudo-R² (total) 0.5245
Fixed Effects
Est. 2.5% 97.5% z val. p
(Intercept) -0.7764 -1.3971 -0.1558 -2.4518 0.0142
s_Cadence 0.9420 0.8870 0.9969 33.5759 0.0000
DayWD -0.2166 -0.4103 -0.0230 -2.1925 0.0283
DayND 1.3800 1.1046 1.6554 9.8210 0.0000
PeriodAfternoon 0.5482 0.3630 0.7334 5.8022 0.0000
Activity4NonWork -0.7817 -0.9181 -0.6452 -11.2296 0.0000
Activity4Driving -0.0684 -0.2238 0.0870 -0.8625 0.3884
s_PA_Percent 0.2893 0.1634 0.4152 4.5035 0.0000
s_PA_Intensity -0.0177 -0.1465 0.1111 -0.2693 0.7877
s_Sleep_Time -0.0928 -0.1676 -0.0179 -2.4287 0.0152
ProfessionSP 0.7743 -0.2048 1.7534 1.5500 0.1211
s_N_MD 0.4721 0.3800 0.5642 10.0420 0.0000
s_N_PD 0.3489 0.2777 0.4201 9.6054 0.0000
s_N_P 0.0698 0.0073 0.1323 2.1900 0.0285
s_Agreeableness 0.2090 -0.2864 0.7044 0.8268 0.4084
s_Extraversion -0.5591 -1.1028 -0.0154 -2.0155 0.0439
s_Neuroticism -0.1907 -0.7244 0.3430 -0.7004 0.4837
s_Openness -0.3048 -0.8123 0.2027 -1.1771 0.2392
DayWD:PeriodAfternoon -0.3690 -0.5963 -0.1417 -3.1814 0.0015
DayND:PeriodAfternoon -0.6703 -0.9527 -0.3878 -4.6506 0.0000
DayWD:s_PA_Percent -0.4303 -0.6003 -0.2603 -4.9620 0.0000
DayND:s_PA_Percent 0.1514 -0.0214 0.3241 1.7177 0.0859
DayWD:s_PA_Intensity -0.0224 -0.1921 0.1473 -0.2586 0.7960
DayND:s_PA_Intensity -0.3932 -0.5699 -0.2165 -4.3613 0.0000
s_PA_Percent:s_PA_Intensity -0.9732 -1.1069 -0.8395 -14.2684 0.0000
DayWD:s_PA_Percent:s_PA_Intensity 0.5549 0.3958 0.7140 6.8342 0.0000
DayND:s_PA_Percent:s_PA_Intensity 0.8440 0.6940 0.9940 11.0275 0.0000
Random Effects
Group Parameter Std. Dev.
Participant (Intercept) 1.0900
Grouping Variables
Group # groups ICC
Participant 24 0.2646

26 Visualization of NHR Stress Model

 

27 Composite Model Figure

##  [1] "T001" "T002" "T003" "T004" "T005" "T006" "T007" "T008" "T009" "T010"
## [11] "T011" "T012" "T014" "T015" "T016" "T017" "T018" "T019" "T020" "T021"
## [21] "T022" "T023" "T024" "T025"